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Dielectric permittivity and Second Harmonic Generation (SHG) studies in the field-cooled mode 
show a linear dependence of dielectric stiffness (inverse dielectric permittivity) on dc bias in PMN-PT 
crystals and SHG intensity in KTaOa :Li at small Li concentrations. We explain this unusual result 
in the framework of a theory of transverse, hydrodynamic-type, instability of local polarization. 



I. INTRODUCTION 

Locally disordered ferroelectrics, such as (1 — 
x)PbMg 1/3 Nb 2/ 303 - xPbTi0 3 (PMN-PT-i), KTa0 3 :Li 
etc., exhibit unusual properties, which find wide 
applications. 2 For example, PMN-PT crystals from the 
morphotropic phase boundary (MPB) compositional 
range possess an extremely high piezoelectric coefficient, 
which is important to transform electric energy to me- 
chanical and back>» 

The high dielectric permittivity, e, and its nonlinearity, 
i.e. strong voltage dependence of permittivity, makes dis- 
ordered ferroelectrics very attractive for applications in 
electrically tunable devices, especially in thin-film form. 
In these and many other applications (e.g. electrooptic 
and electrostriction), disordered ferroelectrics are sub- 
jects to high levels of dc biases (E). In the framework 
of Landau theory, e{E) in cubic crystals is quadratic, at 
small fields, and only even terms are allowed** 

It was shown***i that the giant piezoelectric response 
in BaTiC>3 and PMN-NT is due to "polarization rota- 
tion" , which means that the electric field turns the di- 
rection of polarization and produces mechanical stresses. 
It is important that the potential relief in PMN-PT com- 
positions near the MPB is nearly flat, which allows large 
strains at comparatively small fields. 

These materials have some degree of local disorder 
that is known to produce polar regions of finite sizes. 
The rotation of dipole moments in these regions as well 
as sideway domain wall movements (or the domain wall 
movements in 90° domains) can contribute to dielec- 
tric permittivity in a special way (which we will discuss 
below), especially, when the field-cooled (FC) procedure 
is employed (this helps overcoming barriers within the 
measuring time). 

Another model example considering local dipoles is 
KTa0 3 :Li (KLT). Li in KLT is a dipole centre with 
a large Li-related dipole moment & Electric field aligns 
these dipoles and creates polarization in the direction of 
the field. 9 If one changes the direction of the field in the 
FC mode then the direction of polarization also changes. 



That is, the polarization direction follows the direction 
of the field that should be if the Gibbs statistics works. 
Each of the Li dipoles has some potential barrier, which 
should be overcome that produces viscosity. In ferro- 
electric phases of PMN-PT, some viscosity can arise also 
because of domain wall movements, as well as because 
of random fields. The application of the FC procedure 
allows overcoming this viscosity. 

The ability of an idealized dipole (polarization) to ro- 
tate (without having anisotropy) implies that this dipole 
has a transverse instability. Indeed, the application of 
small fields in the transverse direction results in a huge 
(ideally infinite) response. 10 Random fields, anisotropy, 
and barriers suppress this effect at small external fields 
but some features of the transverse instability can be seen 
at larger fields pi* 

In this paper, we present experiments on locally dis- 
ordered ferroelectrics, PMN-PT. We will show that di- 
electric stiffness (inverse dielectric permittivity), in ferro- 
electric phases of PMN-PT, behaves pseudolinearly with 
the field, above a certain threshold. This behavior does 
not follow the Landau theory for cubic crystals?* The 
suggested in the present paper theoretical description of 
this effect takes into account third-power contributions to 
the free energy of ferroelectric phases. The microscopic 
origin of these contributions is discussed in detail for a 
model solid solution, KLT, and we explain experimental 
data* on the pseudolinear field dependence of the SHG 
intensities in KLT at small Li concentrations. 



II. EXPERIMENT AND RESULTS 

The PMN-PT single crystal samples used in this study 
were transparent plates cut from a flux grown crystals 
prepared at the Physics Research Institute of the Ros- 
tov State Universitypi* Large faces of the samples were 
perpendicular to [100] direction. The sample with com- 
position PMN-PT35 having dimensions of about 2x2x1 
mm 3 was optically polished and annealed at 823 K for 
half an hour in air, in order to minimize residual stresses. 
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FIG. 1: Temperature dependencies of PMN-PT35 dielectric 
permittivity measured in the ZF and FC (E = 0.4, 1 and 3 
kV/cm) modes at different frequences :10 3 , 10 4 , 10 5 Hz (from 
top to bottom). 



Platinum electrodes were sputtered on the large faces. 
Thin (0.05 mm) Pt wires were attached to the electrodes 
by silver paste to connect the sample with the contacts 
of a Linkam HFS91: TS 600 hot-stage used as a sample 
holder. Dielectric measurements were performed at the 2 
K/min heating/cooling rate using a computer-controlled 
impedance analyzer Solartron SI 1260. A blocking circuit 
protected the impedance analyzer for experiments under 
high dc bias. Some details of dielectric studies of other 
PMN-PT crystals have been described elsewhere ,A&i4 
Though we focus mainly on the results obtained in the 
zero field cooling (ZF) and field cooling (FC) modes, field 
heating (FH), and zero field heating after field cooling 
protocols were used as well. 

The temperature dependences of the complex dielec- 
tric permittivity e(T) measured in the ZF mode for 
both PMN and PMN-PT30 single crystals have one dif- 
fused frequency-dependent maximum at about 260 K 
and 400 K ( at 1 kHz), respectively (see for details 
Refs. fiii 4 ]). In order to characterize the PMN-PT35 
sample, we show the e(T) curves measured in the ZF 
mode at different frequencies (Fig.l). Two anomalies 
on the e(T) curve are in good agreement with those 
found in Refs. fiii^iLiLiSi] : a diffused maximum cor- 
responding to the cubic-to-tetragonal phase transition 
and an inflexion at a lower temperature, T\, correspond- 
ing to the MPB [the transition between the tetragonal 
and rhombohedral (or monoclinic, according to differ- 
ent data)i2ii£ii§ii£] phases. The temperature T m of 
the e(T) maximum (435 K at 1 KHz) agrees well 
with the data obtained for PMN-PT crystals of simi- 
lar compositionSfiSii^ The temperature T\ («365-375 K 
at 1 KHz) is comparable with the data for flux-grown 
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FIG. 2: Dielectric stiffness vs. the dc bias for [100] PMN- 
PT35 single crystals measured at different temperatures un- 
der FC protocol at 10 kHz: Rhombohedral (monoclinic) 
phase, 323 K (1), 339 K (2), 353 K (3); Tetragonal phase, 
383 K (4), 413 K (5), 423 K (6); Cubic phase, 459 K (7). Solid 
lines are guides to the eye. 

crystals with similar T m values, 15 but is substantially 
(60-100 K) higher than T x values, reported for PMN- 
PT crystals with similar T m values grown by Bridgeman 
methodii£*i& Such discrepancy seems to be caused by dif- 
ferences in crystal's preparation technique and its origin 
needs additional studies. 

In ZF mode, a pronounced dielectric dispersion is ob- 
served below T m in PMN-PT35 and T m increases with 
the frequency / of the measuring electric field. The 
T m (/) dependence obeys the Vogel-Fulcher relation. 19 
While the value of the attempt frequency f — 5 x 
10 11 Hz is typical for relaxors^ the Vogel-Fulcher tem- 
perature To = 430 K is only slightly lower than T m . 
Though dielectric dispersion decreases below T , it re- 
mains considerable in the ZF mode (Fig.l). In the FC 
mode, the dispersion in the tetragonal phase diminishes 
but, at low E's, remains practically unchanged, close to 
T m . Such behavior is similar to that observed in disor- 
dered PST and PSN-based ceramics and crystalsia^fliSI 
and seems to be due to the presence of the sponta- 
neous (thermally driven) transition from relaxor to a 
mixed (ferroelectric/relaxor) phase containing both the 
nanoscale polar regions and macroscopic ferroelectric do- 
mains of tetragonal symmetry. Both the T m and T\ ex- 
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FIG. 3: Dielectric stiffness vs the dc bias for [100] PMN- 
PT30 and [100] PMN single crystals measured at different 
temperatures under FC protocol at 1 kHz: 350 K (1), 360 K 
(2), 370 K (3), 200 K (4), 250 K (5). Solid lines are guides to 
the eye. 



hibit thermal hysteresis, in agreement with first order 
type reported for these phase transitionsii§iiLi& The de- 
gree of the s(T) maximum diffusion evaluated using 
approach^ is about 20 K, in good agreement with the 
data for PMN-PT31 crystals and PMN-PT35 ceramics^ 

In the FC mode, the application of even a relatively 
small dc bias leads to substantial changes of e(T). Per- 
mittivity in the tetragonal phase decreases with the bias 
while, in the rhombohedral phase, it increases initially, 
but, at higher biases, begins to decrease (Fig. 1). As 
a result, the inflexion transforms into a maximum at 
Ti « 363 K, in good agreement with the data for the 
[001]-oriented PMN-PT crystals from the 0.31 < x < 
0.35 compositional rangepiZ*i& 

As the dc bias increases, the main e(T) maximum, at 
first, becomes sharper and higher but then lowers, dif- 
fuses and shifts to higher temperatures (Fig. 1). These 
changes are in good agreement with the data for the 
[001]-oriented PMN-PT crystals from the 0.3f < x < 
0.35 compositional rang o 17 i 18 as well as with the data on 
disordered PbSc^Ta^Oa ceramics^ The frequency 
dispersion of e and T m reduces dramatically for bi- 
ases exceeding (0.7-1) kV/cm (Fig. 1). Thus, at low E, 
the PMN-PT35 crystal exhibits a relaxor-like behavior, 



while, at high _B's, its properties are similar to ordinary 
ferroelectrics. The details of the dependence of T m and 
Ti on E will be a subject of a separate publication; here 
we focus on the field dependences of dielectric stiffness. 

The main new result, which we want to discuss in the 
present paper, is the pseudolinear field dependence of 
dielectric stiffness in the ferroelectric phases of PMN- 
PT35 above a certain threshold (Fig. 2). Such behavior 
is not unique. For comparison, Fig. 3 presents dielectric 
stiffness vs. E for [100] PMN-PT30 and PMN crystals. 
It is seen that, the pseudolinear dependence is fulfilled for 
these compositions as well, but above a larger threshold. 



III. THEORY 

A. Hydrodynamic model 

At first, we want to discuss experimental data obtained 
for model compounds, KTaOsTi. This case has been 
studied carefully on the basis of first-principles compu- 
tations (see& and references therein) and experiment £ It 
was obtained that Li's are off-center and have a large 
dipole moment. Li-related dipoles (further, Li-dipoles) 
can occupy 6 possible positions with equal energy. Hence, 
each dipole can rotate (FC protocol is assumed). Let a 
dipole occupy a position along z. The application of field 
SE perpendicular to z rotates the dipole. The change 
of polarization can be expressed as: 



6P ± = X ±SE = nfj.SE/E 



(1) 



where n is the Li concentration, fj, the dipole moment, 
E the magnitude of the field along axis z, and x± t ne 
transverse susceptibility: diverges in this simplified 
approach^ but taking into account the random fields and 
anisotropy of the potential relief stabilizes the system^ 
The change of the polarization magnitude in the field 
can be described by the Landau theory based on the soft 
mode concept, which we will discuss below. The free en- 
ergy necessary for the description of the transverse fluc- 
tuations has the formSi 



5F 



dV 



(2) 



where c is constant. One can obtain the longitudinal 
susceptibility from relations between the transverse and 
longitudinal changes of polarization 9 : 



= S J1 = _L A /rxp ) 2 \ = nkBT 

X " SE 2PdE\ { ±j / 8n((u)cf 2 E^ 
E > (3) 

Notice that this susceptibility is due to transverse fluc- 
tuations of polarization. Normally, it should be added to 
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FIG. 4: The field dependence of the FC SHG intensity in 
KLT 9 . Solid lines are guides to the eye. 



the Langevin susceptibility that will be discussed below. 
The integral of this equation gives 2 ^ 



(P - P y = aE 

P>Po 

E>0 




contributions. The evidence of equation in dielectrics 
was found in experiments on SHGS (in magnetic systems 
it was discussed i n 10 i 25 ). 

The SHG intensity, S, is proportional to the average 
square of polarization: 



(6) 



It gives the linear dependence of the FC SHG intensity 
on E at aE ^> Pq . This finding is in agreement with the 
SHG experiment performed for KLT at Li concentration 
1.6% (KLT 1.6%), in the FC mode (see Fig. 4). 



B. Pseudolinear field dependence of dielectric 
stiffness 



The above consideration is true if only the lattice (soft 
mode) contribution to polarization is comparatively neg- 
ligible. Experiments on KLT 2 ' 1 showed that the lattice 
(soft mode) contribution to SHG intensity starts at much 
larger fields than those used in experiment^ In some 
other cases, the contribution of the soft mode may be 
comparable with the contribution of local dipoles, po- 
lar regions or domain walls. We think that this is the 
case in PMN-PT35. Due to low symmetry of ferroelec- 
tric phases (we are talking about the phases obtained in 
the FC mode at E — > 0), free energy includes an odd 
term: 



F = -aP 2 
2 



-AP 3 + -bP A - EP 
3 4 



(4) 



P> P = -a/A 
The equilibrium condition gives 



(7) 



where a — (nfc^T) 2 /167r 2 ((/x) c) 3 and Pq is the FC po- 
larization magnitude at small E (remnant polarization). 
Random fields modify this equation of stated It was ob- 
tained that equation Q is true for the field magnitude, 
which is larger than the random-field magnitude. It is 
important to say that equation Q corresponds to a cu- 
bic contribution to free energjiiS: 



AF = A (P - P f - EP 
P>Po 



(5) 



where A — a/3. 

Notice, that the same conclusions can be got under as- 
sumption of domain wall contributions to dielectric per- 
mittivity. Experimentally, for dielectrics, some violations 
of Landau expansion over the even terms (valid for the 
paraelectric phase only) were discussed for PLZT by J- 
L. Delist; the author— connected this with domain wall 



aP + AP 2 + bP 3 



E = 



P>Po 
Now one can obtain that 



(8) 



X' 1 =a + 2AP + 3bP 2 
P>Pa 



(9) 



At small fields, P = —a/A+\a\ 1 E that leads to the ob- 
served linear field dependence of inverse susceptibility 28 : 



X = 
E>0 



3bP 2 



2A + 6bP 



E 



a 



(10) 



At large fields, one needs to add even terms, and the de- 
pendence becomes quadratic and, then, E 2 ! 3 as it should 



5 




E [kV/cm] 



FIG. 5: The fit of theory (solid line) to experiment (points) 
for PMN [111] obtained in the FC protocol at T=250 K and 
frequency =100 Hz. 



F(e,T c ) = -aP 2 



pP i + - 1 P e -(E + e)P (11) 



where a — a (T — T c ) and e is random field. From the 
equilibrium condition, one gets: 



aP - (3P 3 + 7P 5 = E 



(12) 



Dielectric permittivity can be found by differentiating 



e' (e, T c , E) = 1 



1 dP 

Eq dE 



= 1 



1 



1 



£ a - 3(3P 2 



5 7 P 4 

In the case e = 0, polarization experiences a jump at 
some field Eo, and dielectric permittivity diverges at this 
field. Due to the random fields and distribution of T c 's, 
this jump diffusesi22iSI This can be described by intro- 
ducing a distribution function /(e,T c ) (see the way how 
to determine the distribution function from experiment 
in Ref. 30): 



be in the Landau theory developed for cubic crystals. 
The interval of the fields, where UlUjl is true, depends on 
the coefficients in expansion Q). One can expect this 
effect close to the phase transition temperature, and at 
large A, which is possible when the transverse spatial 
fluctuations of polarization are large. Both conditions 
are satisfied at MPB. 



C. Influence of random fields on the FC 
dependence of dielectric stiffness in PMN 

The influence of random fields on hydrodynamic fluc- 
tuations of polarization results in the insensitivity of di- 
electric permittivity to external fields at small fields>ii 
Here we will discuss another aspect of this problem. 
Experiment shows that, in the ferroelectric phases of 
PMN-PT35, dielectric stiffness pseudolinearly depends 
on the field in a wide interval of the field. Only at very 
small fields, where, perhaps, the time of measuring is not 
enough, there are some deviations from the pseudolincar 
dependence (Fig. 2). Experimental results obtained for 
pure PMNi^ show a more pronounced influence of the 
random fields on dielectric stiffness. It was found that, 
in the [111] direction, dielectric stiffness of PMN has a 
bump while, in the [001] direction, dielectric stiffness de- 
creases with the field smoothly. We suggest considering 
this difference as a result of anisotropy of free energy. We 
take into account that there is a diffuse first order phase 
transition in the [111] field. 14 

We consider the free energy in [111] direction corre- 
sponding to the first-order phase transition between the 
cubic and rhombohedral phases: 



e av (E) = J e{e,T c ,E)f(e,T c )dedT c (14) 

In order to simplify this complex situation, one can model 
inhomogeneous matter by introducing a Gauss type dis- 
tribution of Eq: 

P(E)=P 1 (E) + 

+P 2 {E) f 6{E - E)e^ E °- El V x °f dE (15) 

where E\ is the position of the dielectric permittivity 
maximum, 9(x) is a step function, xq is the width of the 
distribution function, P\{E) and P2{E) are continuous 
functions. Dielectric permittivity can be got by differen- 
tiating ljT5|l : 

e{E) = ai(E) - 
-a 2 (E)eri[(E-E 1 )/x Q } + 

+P 2 {E)e-^ E - E ^' x ^ (16) 

here erf(cc) is the probability integral, aq and a 2 can be 
expressed over the derivatives of Pi(E) and P 2 (E). 

Fig. 5 shows the result of the fitting of our expres- 
sion to experiment. The distribution function width was 
found at 1.97 kV/cm, and E c = 4.27 kV/cm. For the 
sake of simplicity, we took aq and a 2 as constants, and 
P 2 was proportional to E. The fit shows that the nature 
of the bump in the dependence of dielectric permittiv- 
ity on the [111] field can be the diffuse anomaly due to 
the first order phase transition between the cubic and 
rhombohedral phases. 
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IV. DISCUSSION 

The performed experiment reveals that dielectric stiff- 
ness in the ferroelectric phases of PMN-PT from the 
MPB range behaves linearly with E. By integrating i|l(J|) 
we obtained that the dependence of polarization on E 
is logarithmic. This dependence includes both even and 
odd terms although, in the paraelectric phase only even 
terms are expected. In ferroelectric phases, the odd terms 
are allowed due to extrinsic contributions. 26 Notice that 
the paraelectric phase can be a result of ZF procedure. 
Cooling the same sample at finite E can result in the ap- 
pearance of remnant polarization. Our experiments have 
shown that the pseudolinear field dependence of dielec- 
tric stiffness is well observed at temperatures near the 
MPB (where an intermediate monoclinic phase can ap- 
pear stimulating the polarization rotation^), but, in the 
paraelectric phase, the data are not that definite. Thus 
the important feature for the observation of the lineariza- 
tion of dielectric stiffness vs E is a fiat enough potential 
relief, which is indeed the case in PMN-PT compositions 
from the MPB range. Note, that in the PMN-PT30 crys- 
tal, which is outside of the morphotropic phase boundary 
region and especially in pure PMN, the pseudolinear por- 
tion of the \/e(E) dependence is observed at larger fields 
than in PMN-PT35. This is consistent with theory^ 
which showed that the onset of the unusual power law in 
the hydrodynamic model is at the fields, which are larger 
then the random field magnitude. In that study only 
the cubic term was considered in the free energy. In the 
present study we have added the normal quadratic term 
responsible for the polarization of lattice (soft-mode con- 
tribution). We have shown that these two terms together 
result in the pseudolinear field-dependence of nonlinear 
permittivity. 

Our finding correlates with the recent observation 
of the linear field dependence of dielectric stiffness in 
thin films of PMN where the energy barriers between 
metastable states of polar clusters appear smaller than in 
bulk and ceramics^ Approximately linear dependences 
of dielectric stiffness on dc bias were observed in disor- 
dered PST ceramics near the border of stability of relaxor 
vs ferroelectric phased 

At small Pq, the expression obtained for dilute KLT 



provides the uncommon equation of state, P 2 ~ E, which 
resembles that obtained for some magnets^ This be- 
havior has been evidenced in experiments on SHG in 
KTa0 3 :Li for the Li concentrations 1.6%& (Fig. 4), 
which, in our opinion, reflects the fact of strong transver- 
sal spatial fluctuations of polarization. If the ZF proce- 
dure were employed then one would not expect such an 
effect at low T due to the nonergodicity in this system 
arising because of potential barriers (viscosity) . In such 
experiments only the ordinary, Langevin type, response 
of the dipoles should be seen. However, if, due to different 
reasons, a tilt of the Li dipoles is prevented by internal 
random electric fields or stresses then the fluctuations 
become weak and the corresponding contribution to the 
susceptibility vanishes. For example, for the concentra- 
tions larger than 2%, the SHG intensity behaves with the 
field exponentially, ~ exp [(fiE — w) /&bT] (Fig. 4, no- 
tice the logarithmic scale). It can be connected with the 
fact that the Li dipoles are merged into large clusters or 
even domains at these concentrations. A strong Li - Li 
interaction prevents their reorientation. It implies that 
there is energy, w, which is necessary to reorient dipoles. 
The electric field applied decreases this energy linearly. 
Such dependence can be seen even for low Li concentra- 
tions if the sample is strained. There can be also some 
differences due to different protocol in the FC procedure 
(for instance, due to different annealing temperatures). 

Finally, we have observed a pseudolinear field depen- 
dence of dielectric stiffness in [100] PMN, PMN-PT35 
and PMN-PT30 single crystals in the FC mode. We 
have explained this dependence by polarization fluctu- 
ations (rotations) and domain wall movements. 26 Within 
the same approach we managed to explain the experi- 
mental data™ on the pseudolinear dependence of the SHG 
intensities on E in KLT at small Li concentrations. 
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